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Abstract 

A comprehensive microscopic study of the properties of bulk matter at den- 
sities just below nuclear saturation p s = 2.5 ~ 10 14 g cm~ 3 , zero and finite 
temperature and high neutron fraction, is outlined, and preliminary results pre- 
sented. Such matter is expected to exist in the inner crust of neutron stars and 
during the core collapse of massive stars with M > 8Mq 



1 Introduction 

Understanding the phase transition from inhomogeneous to uniform nuclear 
matter is important in the study of a number of astrophysical phenomena. In 
the region of density 0.05p s < p < 0.16,o s , temperature < T < lOMeV and 
proton fraction 0.01 < y p < 0.3, which is expected at the bottom of neutron 
star crusts and during core collapse supernovae, heavy nuclei immersed in a 
fluid of neutrons are expected to become distorted into a series of exotic struc- 
tures known as 'nuclear pasta' pQ in order to minimize the sum of their surface 
tension energy and Coulomb repulsion with adjacent nuclei. These structures 
may extend over distances many times 'normal' nuclear radii w lOfm. 

Such matter could have a significant impact on the dynamics of core collapse 
supernovae. It exists in the outer regions of the collapsed core where neutrinos 
are expected to be trapped. The nuclear pasta phases may affect the neutrino 
opacity of the matter as excitation of collective modes in the pasta offers another 
channel for the transfer of energy from neutrino flux to the nuclear medium [2] . 

In the inner crust of neutron stars the pasta phases are expected to be a 
form of 'soft' condensed matter, that is, in a state between liquid and solid [3]. 
In addition, the external neutron gas is expected to be in a superfluid state, and 
as such its bulk flow will be quantized into vortices, whose cores are comparable 
to the extent of the pasta shapes in size [4] . 

Our theoretical studies are required for the explanation of observations. 
Thermal emission has been observed from seven young neutron stars [5]; to 
constrain models of the cooling of the neutron stars' cores, we must know the 
relevant mechanisms of heat transport in the crust and their efficiency. Quasi- 
periodic oscillations have recently been observed in the tails of flares from soft 
gamma-ray repeaters [B] and have been interpreted as being produced by seismic 
waves in the crust as it relaxes after the significant readjustment that triggered 
the flare. Glitches in the spin down rates of pulsars tell us something about the 
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coupling of the crust to the core [7]. Later in a neutron stars' life, accretion may 
lead to the formation of 'mountains' on the neutron star surface, leading to a 
quadrupole moment and gravitational wave emission: whether the gravitational 
waves can be detected on Earth depends on how big a mountain the crust can 
support [8 . Accretion can also reheat the crust, melting it in layers [9] and will 
also compress the innermost layers of the crust. 

2 Motivation for a New Study of the Pasta 
Phases 

In order to study the structure of the inner crust in an unbiased way, one must 
perform the calculations self-consistently in three dimensions, making no as- 
sumptions about the shape of the nuclear distribution or lattice type (i.e. free 
of the Wigner-Seitz approximation) and no distinction between the nuclear clus- 
ters and external neutron gas. Because of the computationally intense nature 
of such calculations, they have only recently been attempted, by two methods 
in particular; the semi-classical Quantum Molecular Dynamics (QMD) [TT] and 
the fully quantum mechanical Hartree-Fock (HF) method [T^] . QMD allows for 
a large volume (of order lOOfm) to be simulated, and thus longer range (lower 
energy) effects explored. HF is more intensive computationally, so the simula- 
tion volumes are smaller (of order 20-40fm), but effects arising from the shell 
structure of the nuclear clusters and the unbound nucleons which scatter off 
them are automatically included. Both the above studies also neglected the 
band structure of the unbound nucleons, for which one requires general Bloch 
boundary conditions 

We would like to extend the Hartree-Fock study above to a far more com- 
prehensive range of density and proton fraction space, as well as extending it 
to finite temperature thus exploring the properties of matter in core collapse 
supernovae. In doing this we would survey the effects of shell structure on the 
energy and ordering of the shapes, paying particular attention to separating 
physical effects from numerical artifacts, and implementing the general Bloch 
boundary conditions so we can examine the band structure of the neutron fluid. 

3 3d Hartree-Fock Simulation of Bulk, Inhomo- 
geneous Nuclear Matter 

The Hartree-Fock approximation assumes a definite set of quasi-particle states 
with energies ei, ...,6a that are occupied with a probability given by the Fermi- 
Dirac distribution 



where q = p,n labels the isospin states, i the single particle states and /i 
is the chemical potential. Physically we are making the assumption that the 
nucleons in a nucleus or a nuclear configuration move independently of each 
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Free Energy Density vs Baryon Number 
for T = 2.5MeV, y =0,3, n b = 0.083 fm 
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Figure 1: Right: The total free energy density, at temperature, density and 
proton fraction as given, as a function of number of nucleons in the cell (or, 
equivalently, the size of the cell). Left: some nuclear configurations obtained: 
top left - A=460, top right - A=2200, bottom left - A=280, bottom right - 
A=1300. 



other in an average potential created by all the other nucleons. The variational 
principle <5£ , [<f>] = 6($>\H |$) = is used to obtain our approximation to the 
ground state given a Hamiltonian containing a two-body nuclear interaction. 
We use the Skyrme interaction, a zero-range effective nuclear force particularly 
suited to Hartree-Fock calculations. Carrying out the variational procedure with 
respect to the single particle wavefunctions of a single Slater determinant, the 
two body Skyrme potential becomes a one body density-dependent potential 
and the A-body Schrodinger equation for the A-body wavefunction becomes A 
single body Schrodinger equations for the single particle wavefunctions: 



^HF<Pi,q 



J i.q c i,gY-%g 



(2) 



Here, w q is the spin-orbit potential (which we currently set to zero), u q is 
the single particle potential, and m* is the effective mass. See [10] for the form 
of the Skyrme interaction and derived one body potentials. 

We also impose a constraint that the neutron quadrupole moment of the 
nuclear configuration be given by an input value (Q) p - If we didn't do this, we 
would have no control over which of the local minima in deformation space the 
simulation would fall into, destroying self-consistency. 

We make the assumption that at a given temperature and density the matter 
is arranged locally in a periodic structure with an indentifiable unit cell. We 
then take our computational volume to be that unit cell. We impose Bloch 
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boundary conditions: 

^, q (r + f)=e: R - f ^ q (r) (3) 

where T is the translation from the position r to the equivalent positions in 
the adjacent cells, and K is the Bloch momentum covector. We must perform 
one simulation for each value of K within the first Brillouin zone. 

Each unit cell will contain a certain number of neutrons N and protons Z, 
making a total baryon number of A — N + Z. We have freedom to increase the 
cell volume V and number of nucleons A and still describe the same density. 

Thus at a constant density, temperature and proton fraction we must scan 
across different cell sizes, neutron quadrupole moments and Bloch covectors. 
In order to reduce the numerical work, we will actually scan across cell size 
and quadrupole moment with simple periodic boundary conditions {K = 0), 
select the configuration that gives the minimum energy, and then perform the 
calculation for that configuration with K ^ 0. 

We are in the process of performing the calculations outlined above on the 
Cray XT3 super computer Jaguar at Oak Ridge National Lab. A sample of 
preliminary results is given in figure ([1]). 
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